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We present ab initio calculations of the evolution of anisotropic magnetoresistance (AMR) in Ni 
nanocontacts from the ballistic to the tunnel regime. We find an extraordinary enhancement of 
AMR, compared to bulk, in two scenarios. In systems without localized states, like chemically pure 
break junctions, large AMR only occurs if the orbital polarization of the current is large, regardless 
of the anisotropy of the density of states. In systems that display localized states close to the Fermi 
energy, like a single electron transistor with ferromagnetic electrodes, large AMR is related to the 
variation of the Fermi energy as a function of the magnetization direction. 

PACS numbers: 



I. INTRODUCTION 

The free energy of mono-domain ferromagnetic parti- 
cles depends on the relative orientation of the magne- 
tization with respect to the crystal lattice. This mag- 
netic anisotropy results from the combination of Coulomb 
repulsion favoring spin polarization, spin-orbit coupling 
(SOC), and the crystal field breaking the orbital rotation 
invariance. As a result, the orbital moment of magnetic 
atoms and their magnetic anisotropy energy (MAE) de- 
pend strongly on their atomic coordination^. 

The transport counterpart of MAE is anisotropic 
magneto- resistance (AMR), i.e. the dependence of the 
resistance on the angle between the magnetization and 
the current flow. Whereas AMR in bulk was known back 
in the XIX century and is a rather small effect, the re- 
cent observation of AMR in a variety of low dimensional 
systemaS^iSiSiZiSiS i 10 1 1 1 1 1 2 , largely exceeding bulk values, 
has opened a new research venue in the field of spin- 
polarized quantum transport. Very large AMR has been 
reported in planar tunnel junctions (TAMR) with a vari- 
ety of electrode and barrier materials^^&i^. Enhanced 
AMR has also been observed in atomic sized contacts, 
both in the tunnel regime (TAMR) and in the contact 
(or ballistic^) regime (BAMR)H, for Py£, Fe±£, Niii, 
and Coi2. Additionally, GaMnAs islands in the Coulomb 
Blockade regime show electrically tunable AMR (CB- 
AMR) 15 . Thus, a wide range of nanostructures made 
from different materials display enhanced AMR. 

Here wc focus on AMR in atomic-size conductors 
for several reasons. On the one hand, conductance of 
atomic-sized contacts probes the electronic structure of 
the apex atoms. These have coordination different from 
bulk and thus present different orbital and spin mag- 
netic moment^, and enhanced magnetic anisotrop y b 2 i 17 
which might be probed by BAMR. On the other hand, 
nanocontacts allow to study AMR going from the con- 
tact (BAMR) to the tunnel (TAMR) regime in the same 
system, as shown in the case of both Ni and Pyiil&. Ni 
nanocontacts have also been used as electrodes to explore 
the Coulomb Blockade and the Kondo regimes^. 

The crux of the matter is to identify the necessary 



and sufficient conditions to expect large values of AMR 
in quantum transport. Here we consider two different 
transport regimes, coherent and sequential. In the co- 
herent regime we use the Landauer formalism that, at 
zero temperature, relates the zero-bias conductance G to 
the quantum mechanical transmission of the electrons at 
the Fermi energy, G = ^-T(ep). This approach accounts 
for AMR both in the tunneling regime (TAMR)2°i and 
in the contact or ballistic regime (BAMR)^ in the ab- 
sence of sharp resonances near the Fermi energy. In the 
scattering-free case of perfect ID chains, T(ep) is sim- 
ply given by the number of bands at the Fermi energy 
N{cf)- Because of the SOC, N(<£f) for ferromagnetic 
ID transition metal chain o 10 ' 12 : 14 : 17 depends on the an- 
gle 9 between the chain axis and the magnetization, and 
this leads naturally to stepwise G(9) curves. 

However, the idealized scattering- free picture fails to 
account for the experimental results of conductance in 
metallic nanocontacts, for which scattering channels are 
not perfectly conducting 2 ^. According to the scattering- 
free theory, the conductance of atomic contacts of Ni, 
in units of e 2 /h, would be 6 or 7 for Nii^, in quanti- 
tative disagreement with the measured^ 2 - conductance of 
Ni nanocontacts around 3e 2 /h. The same applies to Fe, 
Co and Pt. Scattering definitely affects d-bands, which 
suffer the so called orbital blocking 2 ^. 

Here we present calculations of BAMR plus scattering. 
This approach also permits to calculate the crossover 
from BAMR to TAMR. Wc find that in the coherent 
regime large AMR is related to the orbital polarization 
of the current. TAMR has been linked to the anisotropy 
of the density of states at e^, which turns out to be large 
in Ni chains. Unexpectedly, this does not lead to a large 
value of TAMR, the reason being that the current is not 
orbitally polarized in this limit. 

In the sequential regime, valid to describe systems 
that feature transport through resonant levels of width 
r smaller than the temperature kT—, we find enhanced 
AMR, regardless of the orbital polarization of the cur- 
rent, if the chemical potential cf of the ferromagnetic 
electrode crosses a resonance as it varies due to change 
of the magnetization angle. This situation occurs in a sin- 
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FIG. 1: ID Model of Ni nanocontact: Elastic electron scat- 
tering in the contact region of a real nanocontact is mimicked 
by a gap d > a between two semi-infinite ID Ni chains with 
lattice spacing a. 



gle electron transistor with ferromagnetic electrode a 15 ' 25 . 
Resonances might also occur in the tip atoms of Ni 
nanocontacts in the tunneling regime^. 

This paper is organized as follows: In Sec. |TT]we intro- 
duce the model system and the theoretical method for 
the calculation of AMR in magnetic nanocontacts. In 
Sec. IIIII we calculate the AMR in the contact or ballistic 
regime (BAMR). In Sec. HVlwe treat AMR in the coher- 
ent tunneling regime (TAMR), while in Sec. |V]we treat 
AMR in the sequential tunneling regime. Finally, in Sec. 
I VII we conclude the paper summarizing the main results. 



II. MODEL AND METHODOLOGY 

As a model system for ferromagnetic nanocontacts we 
consider two semi-infinite Ni ID chains with lattice pa- 
rameter a, separated by a gap d > a, as shown in Fig. 
[TJ This model shares most of the relevant features with 
realistic nanocontact models, like e.g. the low coordi- 
nation of the tip atoms of the two electrodes and elas- 
tic electron scattering due to the gap. On the other 
hand the one-dimensionality and the resulting rotational 
invariance of our model simplify considerably the cal- 
culations of the SOC and the interpretation of the re- 
sults. Such one-dimensional models have been employed 
before to study fundamental properties of atomic-size 
nanocontact a 27 ' 28 ' 29 . 

We calculate the electronic structure of the system us- 
ing a combination of density functional theory (DFT) 
in the local spin density approximation (LSDA) and a 
Green's function technique to account for the fact that, 
when d ^ a, the system is not translationally invariant. 
We split the system into 3 regions, left (L) and right (R) 
electrodes, described as semi-infinite Ni chains, and the 
central region (C) containing the 3 innermost atoms of 
each electrode, as shown in Fig. [TJ 

The electronic structure of both the electrodes and 
the central region is described using effective one-body 
Hamiltonians obtained from ab initio calculations, per- 
formed with CRYSTAL0322 on the LSDA level, and 
using a localized atomic orbital minimal basis set. 
CRYSTAL03, which docs not include SOC, yields spin- 
polarized solutions along an arbitrary axis with majority 
and minority electrons. The SOC term Hgo = AZ • S 
is added to the converged self-consistent LSDA Hamilto- 



H = H\ 
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This post self-consistent approac h 14 ' 31 is justified in the 
case of Ni, for which the SOC is much smaller than the 
exchange splittings and the bandwidths. We take A = 
70meV for the Ni 3ci-orbitals. 

The Green's function of the central region is obtained 
by means of the so-called partitioning technique: 



G C (E) = {E-H C - t L {E) - t R (E)Y 



(2) 



where He is the total Hamiltonian (including SOC) of 
the C region, and El and £_r are self-energies that take 
into account the coupling of the central region to the two 
electrode a 32 ' 33 . 

In the coherent regime, i.e. at low temperatures and 
small bias voltages when inelastic scattering events can 
be neglected, we use the Landauer formalism to calculate 
the conductance of the system which is obtained from the 
transmission function. The transmission function in turn 
can be calculated by means of the Caroli expression from 
the Green's function of the central region, eq. (J5J), and 
the so-called coupling matrices, Tl(E) := £(Ez, — St) 
and t R (E) := i(t R - t R ), of the electrodes^: 



T(E) =Tr [G C (E)T L (E) & c (E)f R (E) 



(3) 



The zero-bias conductance is then given by G — %-T(ep). 
The orbital projected density of states p a and the density 
of states of the central region p can be calculated from the 
Green's function: p a (E) = — -lm[G aa (E)] and p(E) = 

-ilmTr[G c (£)]. 



III. CONTACT REGIME 
A. The ideal chain 

Magnetic anisotropy comes from the combined action 
of both, the crystal field that breaks the orbital rotational 
invariance, and the atomic SOC term, that couples the 
spin polarization to the orbital degrees of freedom. The 
electronic structure of the ideal one dimensional Ni chain 
(d = a) presents a number of common features with 3d 
and 4d transition metal a 17 ' 35 and permits to understand 
transport results for d ^ a. The bands close to the Fermi 
energy arc formed by s and d orbitals. In the absence of 
SOC, rotational invariance around the chain axis permits 
to classify the d orbitals according to the projection of 
their angular momentum along the chain direction, m z . 
On top of that, a weak crystal field splits the otherwise 
degenerate d levels into two doublets E\ (linear combi- 
nation of states with m z = ±1) and Ei (m z = ±2) and 
a singlet A\ (m z = 0), which is hybridized with the s 
orbital. The orbital degeneracy of the doublets is kept 
by the bands of the chain, as long as SOC is not present. 
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FIG. 2: (a) Transmission for an ideal infinite Ni chain (d = a, 
upper curves) and for two semi-infinite Ni chains seperated 
by d = 1.3a (lower curves) for magnetization parallel (black) 
and perpendicular (grey) to the chain axis, (b) Zero-bias 
conductance G as a function of d for magnetization parallel 
(black) and perpendicular (grey) to the chain axis, (c) BAMR 
(grey boxes) and OPC (black circles) as a function of d. (d) 
G as a function of 9 for different values of d/a: 1.0 (full black 
boxes), 1.1 (grey circles), 1.2 (black full circles), 1.3 (grey 
triangles) and 1.4 (black full triangles). 



The bandwidth of the Ei is significantly smaller than 
that of E\ due to the smaller overlap of the E2 orbitals. 
As a result, the E 2 bands yield a higher density of states. 

The combined action of SOC and magnetism alters 
this situatio n 14 ' 17 . When the magnetization is pointing 
perpendicular to the chain axis (9 = 90°), SOC acts as 
an effective magnetic field acting over L x that has to 
compete with the L 2 -like terms arising from the crystal 
field, which happens to be dominant. As a result, the 
bands for 9 = 90° look very similar to those without SOC, 
except in the points where bands with m z o~ and m z ± 
1, <7T-1 intersect, which are far away from tp in the case of 
Ni. Therefore, for 9 = 90° the effect of SOC on transport 
is negligible. In contrast, when magnetization is pointing 
along the chain axis (0 = 0°), SOC shifts the bands by an 
amount Xm z a, where m z and a are the projection of the 
spin and orbital momentum along the chain. As a result, 
the E2 and E\ orbital doublets are split so that one of the 
2 minority E2 bands is shifted below the Fermi energy, 
compared to the 9 = 90° case. This can be seen in the 
stepwise curves in Fig. 2a, that correspond to Xji (E) = 
T(E, 6 = 0°) and T ± (E) = T(E,6 = 90°) for the ideal 
chain. At the Fermi energy, T±(E) ^ Tu (E). This change 



is responsible for BAMRM, defined as BAMR= - 
where AG = G± - G\\. 

The interplay between SOC and magnetization results 
in a non-zero orbital moment density along the magne- 
tization direction. The largest orbital moment occurs 
when 9 = 0, i.e. when the magnetization is along the 
chaini 7 -. The orbital polarization current (OPC) defined 



as 



V T 
OPC = - m ' 



T(E) 



(4) 



where T m is the transmission of the d-orbitals with m = 
±2 or m = ±1 along the chain direction, vanishes when 
9 = 90° but is non-zero when 9 = 0. Interestingly, there 
is a perfect one-to-one correspondence between the OPC 
and the BAMR in the case of the ideal chain without 
scattering. It is also apparent that the existence of an 
orbital magnetic moment is not a sufficient condition for 
having a non-zero OPC, very much like spin-polarization 
does not necessarily imply a spin-polarized current*^. 



B. The effect of weak scattering 

Now we see how elastic scattering, controlled with the 
chain separation d, affects BAMR. The stretched bond 
mimics the contact region. This perturbation preserves 
the axial symmetry of the ideal chain but introduces scat- 
tering. As a consequence T(E), shown in Fig. 2a, is not 
quantized anymore, as expecte d 22 ' 23 and yet the BAMR 
(Fig. 2c) is close to that of the ideal case for values of 
d/a < 1.4. Relatedly, the G{9) curve is not stepwise (as 
is the case of the ideal chain) anymore when scattering 
is included (Fig. 2d). On the other hand the G(9) curve 
is also different from bulk behavior where G(9) cx cos 2 9. 
The quantized step in the ideal case (d = a) that cor- 
responds to the critical angle at which the E2 band is 
pushed below the Fermi energy^, becomes progressively 
smoother as the gap between the chains increases. Our 
G{9) curves including scattering agree with those of the 
experiments^. This is one of the important results of 
the model. 

As d increases, the scattering increases and G goes 
down but interestingly, the BAMR signal first increases 
slightly for d/a < 1.3 before finally going down with in- 
creasing scattering. The initial increase in BAMR is re- 
lated to the initially stronger decrease of the contribution 
of the ^li-channel to the conductance. The Ai-channcl 
is not affected by the SOC in contrast to the E2 chan- 
nels that are mainly responsible for the BAMR signal. 
The decrease of the BAMR signal for larger values of 
d is expected within the framework of our model, since 
the relative contribution to the conductance of the d- 
channels compared to the s-channel decreases as the gap 
opens. The reason for this is the shorter spread of the 
(i-orbitals compared to the s-orbitals^ 3 -. Relatedly, the 
OPC (Fig. 2c) also decreases as d increases. Remov- 
ing the contribution of the s-channcl to the conductance 
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FIG. 3: Tunnel regime (d = 4a): (a) Transmission function 
for magnetization angles 8 — 0° (black) and 9 = 90° (grey), 
(b) Zero-bias conductance in dependence of the magnetiza- 
tion angle 9. (c) DOS projected onto tip atom as a function 
of energy for 9 = (black) and 6 = 90° (grey), (d) DOS 
projected onto tip atom at Fermi level as a function of 9. 



would thus enhance BAMR. This could be accomplished 
e.g. by oxidation of the contact 3 ^. 



IV. COHERENT TUNNELING REGIME 

In this section we study the anisotropic magnetore- 
sistance in the regime of weakly coupled semi-infinite 
chains. In Figs. 3a and 3c we plot the Landauer trans- 
mission T(E) (calculated from the Caroli expression eq. 
(J3J)) for d = 4a, definitely in the tunnel regime, and the 
density of states (DOS) projected onto the tip atom of 
a semi-infinite Ni chain both for 8 = 0° (pn(E)) and 
8 = 90° (p±(E)). The very small transmission is domi- 
nated by the s-channel, and therefore quite independent 
of 8. In contrast, the DOS is very different for 8 = 0° 
and 8 = 90°. The two peak structure around ep for 
8 = 0° is related to the split E% bands, which merge 
when 8 = 90°. In Figs. 3b and 3d we plot the zero- 
bias conductance G(8) and the DOS at the Fermi level 
p(ep, 8) as a function of 8. Whereas the maximal change 
in the conductance is smaller than 1%, the change in the 
DOS exceeds 200%. This challenges the simplistic link 
between DOS and tunnel conductance. 

In the tunneling regime the Landauer formula can be 
rewritten as (see appendix): 

4 e 2 

GTunncl = ~Y. I ^ I Vq M (5) 

a,p 

where V a p is the matrix element of the Hamiltonian con- 
necting the a and (3 atomic orbitals of the tip atoms of 
the two Ni chains and Pa' R ( e F) is the orbital-resolved 
DOS at the Fermi energy, i.e. the DOS projected onto 



an atomic orbital a of a tip atom. Using this expres- 
sion, the conductance calculated in Fig. 3a from the 
Caroli expression is indeed nicely reproduced. Note, that 
the standard approximation by which the conductance is 
proportional to the product of the DOS of the tip atoms, 
G K E Q Pa( e F) Pp ( € f) is obtained only if the V a p 
matrix is assumed to be proportional to the identity, i.e. 
the tunneling matrix elements are assumed to conserve 
the orbital index and to be equal in size. However, this is 
far from being the case when d and s orbitals are involved. 
In fact, in the case considered here the conductance is 
completely dominated by the V s , s term, for which the 
orbitally-resolved DOS p% ,R is essentially independent of 
8. As a result, the strong dependence of the global den- 
sity of states on 8 is not followed, in this case, by a strong 
dependence of the conductance on 8. Notice that since 
the transmission is dominated by the s channel both the 
orbital polarization of the current and the AMR are neg- 
ligible. In general, an anisotropy in the DOS is not a 
sufficient condition to have AMR. 

The small variation of G(8) in Fig. 3b can be traced 
back to the variation of as a function of 8 and the 
non-flat p^' R (E). In Fig. 4a we plot cf(8) for a semi- 
infinite Ni chain. Aep = ej?(0 o ) — £f(90°) can be as large 
as 10 mcV. This change leads naturally to the second 
scenario for enhancement of the AMR, considered in the 
next section: in a situation of resonant transport through 
the change of the chemical potential ep as a function of 
the magnetization direction can result in a large variation 
of G, regardless of the degree of orbital polarization of 
the current. It has been recently suggested that these 
resonances could arise as localized tip states in Ni wires 
thicker than those considered here^. 



V. SEQUENTIAL TUNNELING REGIME 

In this section we consider a different scenario, moti- 
vated by recent experiments^ and by the remarks at the 
end of the previous section. We study a single electron 
transistor (SET) with Ni electrode o 19 i 37 ' 38 and a non- 
magnetic central island (CI) with a discrete electronic 
spectrum. The CI is weakly coupled to the electrodes, so 
that the levels acquire a broadening T. The position of 
these levels can be electrically tuned with a gate. When- 
ever a level of the CI is in resonance with the Fermi en- 
ergy the zero bias conductance of the system has a max- 
imum. We assume that both the level spacing of the CI 
states, AE, and the charging energy Eq are much larger 
than the temperature fc^T which is larger than T. Under 
these conditions, the system is in the Coulomb Blockade 
regime. 

In equilibrium the chemical potential of the central 
island and that of the electrodes must be the same^i: 
ep(8) = Ec + ejv + eVc where N is the number of elec- 
trons in the CI that satisfies this condition and is 
the energy level occupied by the last electron. ^From 
this equation we immediately see that the charge state 
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FIG. 4: (a) Variation of the Fermi energy of the Ni chain as 
a function of 8. (b) G(Vg, 8) for a SET device coupled to the 
Ni chains. The curves are vertically shifted, (c) G(5meV, 8) 
in logarithmic scale. 



of the central island can be controlled both with the 
gate and with the orientation of the magnetization of the 
electrode a 15 ' 25 . This effect is reminiscent of the so called 
magneto Coulomb effect, in which the chemical potential 
of the electrode is varied with the intensity of the ap- 
plied field 3 ^. Here the chemical potential is changed by 
rotating the applied field. 

In the Eq > kT > T situation, the linear conductance 
of the SET can be obtained using either the finite temper- 
ature Landauer approach 3 - 3 - or the sequential transport 
theory 2 ^: 

G =^cosh- 2 (A/2 fcB T), (6) 



where A = E N {V G ) + - e F {6). In Fig. 4b we 
plot G(V G ,8) for a SET with Ni electrodes. We take 
kT = 5T = 0.5 mcV. The gate is chosen so that, for 
9 = the conductance is maximal. As 9 changes the 
chemical potential of the electrodes moves away from the 
peak. In Fig. 4c we plot G{9) for Vq corresponding to 
the vertical line in Fig. 4b. Notice the logarithmic scale 
and the huge AMR, which might have practical applica- 
tions. Notice that crossing the conductance peak, either 
by gate application or magnetization rotation, implies 
charging the CI by one electro n 15 ! 25 . The results of Fig. 
4b assume that T is independent of 9, which is true as 
long as the resonant level is not coupled to the E2 and E\ 
bands. The height of the G(Vg) curves would depend on 
9 otherwise. In principle, a complete characterization of 
the G(Vc, 9) curve would yield the £f(6) and T(9) func- 
tions, which would provide valuable information of the 
electronic structure of the electrodes. 



VI. SUMMARY AND CONCLUSIONS 

We have presented ab initio quantum transport calcu- 
lations of Ni nanocontacts as a function of the magneti- 
zation direction, 9, going from the ballistic to the tunnel 
regime. We have shown that AMR is unrelated from 
quantization of conductance, which is an artifact of the 
scattering free calculations and not expected in transition 
metal nanocontacts. We also show that a large variation 
of the density of states at €f as a function of 9 is not 
a sufficient condition for large AMR. We identify two 
sufficient conditions to obtain largely enhanced AMR in 
quantum transport. First, in the coherent regime (con- 
tact and tunneling), large AMR is related to a large de- 
gree of orbital polarization of the current, for a selected 
direction of the magnetization. Second, in systems with 
resonances close to €f , as it happens in single electron 
transistors with ferromagnetic electrodes, large AMR is 
related to a large variation of the chemical potential cf 
of the electrode as a function of 9. We report an ab initio 
calculation for this quantity. These findings shed light on 
the choice of materials and the design of nanostructures 
with enhanced anisotropic magnetoresistance. 

We acknowledge R. Aguado, L. Brey, E. Tsymbal, E. 
Tosatti and C. Untiedt for useful discussions. We ac- 
knowledge Spanish MEC and Generalitat Valenciana for 
funding grants MAT2007-65487, Ramon y Cajal Pro- 
gram, GV-ACCOMP07/054 and Consolider CSD2007- 
0010. 



APPENDIX A: DERIVATION OF TUNNELING 
FORMULA 

For completeness, we derive eq. ([5]) from the Landauer 
formalism in the limit of weak coupling between the elec- 
trodes. Eq. can also be obtained from Kubo formula 
(see e.g. the book by Mahan 4 ^ Sec. 9.3). Derivations 
similar to ours can be found in the literatur e) 33 ' 41 . 

We consider two semi-infinite electrodes L and R with 
atomically sharp tips separated by a distance d, as shown 
in Fig. Q] We label the tip atoms of the left and right lead 
and 1, respectively Now the Green's function projected 
onto tip atom is given by: 

Go(E) = (E-H Q -£ L (E)-£ R (E))-\ (Al) 

where S_l is the self-energy representing the rest of the 
left electrode without tip atom while presents the 
self-energy of the entire right electrode including the tip 
atom 1. Thus the right self-energy can be expressed by 
the Green's function of the isolated right electrode 
and the coupling V between the left and the right tip 
atom as: 

£ R = Vg?V^ (A2) 

In the tunneling regime, i.e. for d >> a, when the cou- 
pling V becomes very weak, the contribution of the right 
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sclf-cncrgy to Go can be neglected, so that Go becomes 
equal to the Green's function of the isolated left lead pro- 
jected onto the tip atom, : 

6 (E) ^l(E))- 1 = gk(E). (A3) 

The Caroli expression^! for the Landaucr transmission 
through the tip atom thus becomes: 

T(E) « TT[g^E)t L (E)(g^(E)t R (E)}. (A4) 

The coupling matrix of the right lead Tr can be re- 
written in terms of the spectral function of the iso- 
lated right lead projected onto the tip atom, af' := 

f fl : :=i{± R -t R ) =Vaf Vl (A5) 

The first three terms in eq. (|A4[) are computed using 
the algebraic identity gftt L (gffl = i(g£ - (ffo ) f ) = a o 
where is the spectral function of the isolated left lead 
projected onto the tip atom 0, we find for the transmis- 
sion in the tunneling regime: 

T(E) « Tr [a L G {E) T R {E)] = Tr [a L {E) V a? (E) V\ 

(A6) 



Thus the zero-bias conductance which is given by the 
transmission function at the Fermi energy, can be ap- 
proximated in the tunneling regime by: 

2 2 
G = jxT(e F )^^-xTr[a^(e F )Vaf(e F )V^] 

e 2 x 

= T X, aZ a >(e F )V a ?{ia%fy(eF)VZfi>, (A7) 

ct,a',0,0' 

where in the last step we have labeled states on the left 
tip by a and a' and states on the right tip by (3 and 
(3' . The spectral functions a L and a fl are diagonal in the 
basis of eigenstates of the isolated left and right lead, and 
the diagonal elements yield the DOS projected onto the 
eigenstates: a^ a = 2pa and a§p = 2p^ where a and (3 
now label the projections of the eigenstates onto the tip 
atoms. Thus we obtain eq. ([5]): 

G^^E P L a {e F ) V aP pf(e F ) V* . (A8) 

a,/3 

Notice that this result relates the tunnel conductance to 
the product of the orbital-resolved DOS of the electrodes, 
as opposed to the total DOS. 
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